Angular dependence of the radiation power of a Josephson STAR-emitter 
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We calculate the angular dependence of the power of stimulated terahertz amplified radiation 
(STAR) emitted from a dc voltage applied across a stack of intrinsic Josephson junctions. During 
coherent emission, we assume a spatially uniform ac Josephson current density in the stack acts as a 
surface electric current density antenna source, and the cavity features of the stack are contained in 
a magnetic surface current density source. A superconducting substrate acts as a perfect magnetic 
conductor with H\\^ac = on its surface. The combined results agree very well with recent exper- 
imental observations. Existing Bi2Sr2CaCu208+5 crystals atop perfect electric conductors could 
have Josephson STAR-emitter power in excess of 5 mW, acceptable for many device applications. 

PACS numbers: 07.57.Hm, 74.50.+r, 85.25.Cp 



At present, broad-band terahertz (THz) electromag- 
netic (EM) waves generated from femtosecond laser 
pulses and monochromatic THz waves generated by laser 
mixing, parametric resonance techniques, and quantum 
cascade lasers, etc., are the most common THz sources [ij. 
But these techniques are not cost effective in the "THz 
gap" region 0.1-10 THz required for many important ap- 
plications. After many years of effort, by application of a 
static voltage across the intrinsic Josephson junctions in 
Bi2Sr2CaCu208+5 (BSCCO) mesas, coherent radiation 
power of 0.5 fiSSf was achieved [2]. The same technique 
on different samples since led to radiation power of 5 /iW 
and an output efficiency of 3 x 10~^[3i] . This is an exciting 
development, as the "THz gap" may soon be filled. 

In both experiments, the onset of the intense, coher- 
ent THz radiation occurs in or near to the region in the 
current- voltage (/ — V) characteristic of negative differ- 
ential resistance (NDR), as for the Gunn effect in GeQ. 
In the absence of more precise information, we assume 
that the stack of Josephson junctions acts partly as a 
cavity, in order to amplify the coherent radiation at the 
fundamental Josephson angular frequency uoj and possi- 
bly its harmonics, and partly as a conductor with an ac 
Josephson current, Jac{x' ^ t) = J{x') Yl^^i sm{nujt)^ 
where ujj = 2eVo /{Nh) is the Josephson angular fre- 
quency, Vb is the dc voltage applied across the coherent 
stack of N ^ 10^ junctions, e is the electric charge, h 
is the Planck constant divided by 27r, and is the rel- 
ative amplitude of the nth harmonic of the intrinsic ac 
Josephson current. The effective radiation sources at the 
mesa edges are respectively the surface magnetic current 
Ms arising from the electric field in the cavity gener- 
ated by the non-uniform part of Jad^' ^t)[6]^ and the 
surface electric current Js arising from the uniform part 
of Jac{x',t)^ as sketched in Fig. 1 (a). These sources are 
obtained from the Faraday and Ampere boundary condi- 
tions, respectively [11,0], as used previously 0,0,0, [HI- 
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FIG. 1: (a) Sketch of a mesa with Js and Ms surface current 
sources, (b) Mesa during coherent emission with applied dc 
/o, Vo, and lac confined to it. Curve c is the integration path 
for the Ampere law boundary condition. See text. 



We first assume the mesa is suspended in vacuum. In 
Lorentz gauge, the vector potentials from the respective 
radiation sources are [1,0, 



A{x,t) = 



47r 



R 



(1) 



R = \x — x'\^ k'^ = nujj^/elc = kn\f^^ inside and out- 
side the mesa, e is the dielectric constant of the mesa, 
F{x^t) is given by A{x^t) with eo and Jsni^') replaced 
with /io and Msn{x^)^ respectively, and eo, /io, and c are 
the vacuum dielectric constant, magnetic permeability, 
and speed of light, respectively; Jsn{'^') and Msn{'^') 
are the electric and magnetic surface current densities 
corresponding to the nth harmonic of Jad^'^t). 

We first consider cylindrical mesas of radius a ~ 50 fim 
and height h ^ 1 /im. In cylindrical {p'^(j)'^z') coordi- 
nates, the sources Jsn{x') and Msni^') are 



Jsn{x') = z'^^^ar]j{z')8{p' - a), 



(2) 



Msn{x') = ct)'Sn,i—acos{(l)' - (j)o)r]M{z)S{p - a),{3) 
where Eq = EoJi{k[a) is the effective electric field ampli- 
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tude at the mesa edge, k[a = 1.8412 satisfies J[{k[a) = 
in order that H(f){p' = a) = for the lowest trans- 
verse magnetic (TM) cyhndrical cavity mode, TMf^^gQ, 
Jm{x) is a standard Bessel function, and for no sub- 
strate, rij{z') = r]M{z') = T]{z') = Q{z')Q{h — z')^ where 
Q{x) is the Heaviside step function. The higher en- 
ergy cavity mode energies are not integral multiples of 
one another 0. Assuming the fundamental frequency of 
a non-uniform Jac{x',t) excites a cavity mode, higher 
harmonics of it will not, so we take Msn ^ ^n,!- 
Since h << a^r is the smallest length in the problem, 
r]{z^) hS{z') for no substrate. Outside the mesa, 
we write A{x^t) and F{x^t) in spherical coordinates 
(x, z) = r (sin cos ^, sin 6 sin ^, cos 6) . After averaging 
over ^0 and in the radiation zone r >> a far from the 
mesa, the emitted power per unit solid angle is 



dP 



r/a^oo 



Zo{Jjvki 
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\nanSi{d)Jo{nkg) 



+a(0) (cos^ OJlikg) + Jlikg) 



where a{e) = i|^o5i^(^)lV(^o Jj)', J±{z) = [J2{z) ± 
Jo{z)]/2^ kg = kiasinO^ Zq = \/eo/ /hq is the vacuum 
impedance, (6) = S:^{6) = 1 for no substrate, and the 
details are presented elsewhere [12]. In Fig. 2, plots for 
no substrate of the intensity I{0) ex dP{0)/dVt versus 
at the fundamental with kia = 1.8412/y^ from separate 
Msi and Js sources are shown by curves (A) and (B), 
respectively, and the combined fundamental output with 
a{0) = 0.6 is shown by curve (C). I{0) at the second ac 
Josephson current harmonic with k2 = 2ki is shown by 
curve (E). 

We now consider a rectangular mesa of width length 
^, and height h with no substrate. To fit experiment, we 
assume TM^qq modes, which oscillate in position with 
integral multiples of half-wavelengths along the mesa 
widths [2, 0, 0, [l^. In rectangular source coordinates 

Jsn{x') 



-Vj 



cr=d= 



Msn{x) = ^^r]M{z)sm[n{x - Xn)7r/w] 



X ^ cr[y'fa{x', y') - x'g^{x', y')], 

cr=± 

w5{x' ^ (jw/2)Q[{l/2f - {y')\ 
l5{y'^cjl/2)Q[{w/2f-{x')\ 



(6) 
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where the TM^qo cavity mode energy is degenerate for 
—w/n < Xn < w/n. 

We treat Msn in two models. In Model I, the 
average output power is taken to be {P{xn))i = 
l[P{0) + P{w/n)][6]. In Model II, we let {P{xn))ii = 
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FIG. 2: (color online) Polar plots of the radiation zone inten- 
sities 1(0) from Eq. (|4| in arbitrary units versus (degrees) 
of magnetic and electric source fundamental outputs [(A),(B)] 
and the combined output with a(0) — 0.6 with no [(C), (E)] 
and superconducting [(D), (F)] substrates, of the fundamen- 
tal with k'la = 1.8412 and the second ac Josephson current 
harmonic, respectively, with k2 — 2ki for a cylindrical mesa. 



J^^^^dxnP{xn)n/2w. Then, in spherical coordinates, 
the time-averaged power per unit solid angle in the radi- 
ation zone is 



dP _ Zpjjjvkiy 

dn ^/^^^ 1287r2 
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Xn = {knW / 2) sin COS (j)^ Yn = {kni / 2) sin 



I, II, an{0) = \EonS}fmy{2ZoJjf 



sin(/), where 
, V = wih^ 
and El^iO^cj)) are given 



and the Ci(<9, 
else where 1 12]. 

Plots of I (0^0) oc dP{0^0)/dQ in arbitrary units ver- 
sus 6 in degrees at (/) = from Eq. (|9|) with S:^{6) = 
S^{e) = 1, £ = 20^/3, k[ = tt/w = ki^e, k2 = 2ki 
are given by curves (A) and (B) in Fig. 3[l2[. In 
experiment [3|, the maximum fundamental /(^, 0) is gen- 
erally at (9max ^ 30 - 40°, and 7(90°, 0) = 0. Although 
curve (A), obtained from Msi alone, yields 7(90°, 0) = 0, 
it has 6>max = 0. Hence, it is necessary to add the Jsi 
source. Curve (B) is obtained for a(0) = 0.20, 0.40 for 
i = 1 and II, respectively. It has 6>max ^ 40°, nearly as ob- 
served, but also yields a large 7(90°, 0) value, unlike the 
observations. The corresponding second harmonic inten- 
sity is shown in curve (C), which also shows 7(90°, 0) 7^ 0, 
unlike present experiments Q. This unexpected combi- 
nation of 7(90°, 0) = and ^max ^ 30° led us to con- 
sider the effects of existing BSCCO substrates [H, As 
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FIG. 3: (color online) Polar plots of the radiation zone inten- 
sities /(^, 0) from Eq. (|9)) in arbitrary units versus (degrees) 
with no [(A) Ms source only, (B) combined, (C) n = 2, Model 
II] and superconducting [(D), (E), n = 2, Model I] substrates, 
of the fundamental with k[w — tv and the second harmonic, 
respectively, with = 27r normal to the length (0 = 0) of 
a rectangular mesa with i = 20w/3. 

sketched in Fig. 1 (b), during coherent Josephson radia- 
tion, the ac Josephson current is essentially confined to 
the mesa. With only a dc surface current density cx Iq 
and 5|| cic(^) = beyond the skin depth (^ 0.15 jam 
) inside the BSCCO substrate due to the ac Meissner 
effect, the Ampere-Maxwell boundary condition forces 
H\\^ac{t) = just above the BSCCO substrate [5]. This 
corresponds to a perfect magnetic conductor (PMC) sub- 
strate, with the effective image sources opposite to those 
of a PEC substrate^, J^]. Thus, for a BSCCO substrate, 
we replace r]j{z') and r]M{z') in Eqs. ([2]) and ([5]) by 

n.{z') = ^{z')-^{-z')=sgn{z')Q[h^-{z'n{ll) 

In the radiation zone, h « a^r and h << w^i^r for both 
mesa types, so we assume <C l//cn for the relevant n. 
Expanding e^^"^-^ / R in Eq. ([T]) for small z' ^ 

Si{e) = S^{0) ^Zo. -iknhcose 6(90° - e). (12) 

Plots of I{0) in the radiation zone for a cylindrical 
mesa with kia = 1.8412/^^ and k2 = 2ki atop a super- 
conducting substrate are shown by curves (D) and (F) 
in Fig. 2. Near 90°, a superconducting substrate has a 
drastic effect on the power emitted from the Js source. 
Plots of /(^, 0) for a rectangular mesa on a superconduct- 
ing substrate in the radiation zone of the fundamental in 
both models and of the second harmonic in Model II, 
are shown respectively by curves (D) and (E) in Fig. 3. 
Since intense higher harmonics can arise only from the 
Jsn in cylindrical mesas, the study of cylindrical mesas 
could provide valuable information regarding the primary 
radiation source. 



Preliminary experimental fundamental /(^, 0) results 
for mesas with w = 60 /im, £ = 400 /im, and h = 1 /im 
are in good agreement with curve (D) shown in Fig. 
3^. The effect of a superconducting substrate is cru- 
cial for ~ 90° (nearly parallel to the substrate), as 
shown in Figs. 2 and 3. More importantly, a super- 
conducting substrate has a drastic effect upon the radi- 
ation power. Since ki = n/w for rectangular mesas [s^, 
at the fundamental 6>max, 40° from curve (B) in Fig. 3, 
|5'/'^(6>max)P = [hki cos(40°)]2 ^ 1.6 x 10-3. Hence, re- 
placing the superconducting substrate by an insulating 
one could enhance the power output by 600. Replacing 
it by a standard PEC could further quadruple it(6|, Il2| . 
Since output power of 5 jaW was achieved for rectangular 
mesas on superconducting substrates [3], coherent THz 
radiation power in excess of 5 mW, acceptable for many 
applications, might be attainable from existing BSCCO 
samples. Since the frequency range of the coherent radi- 
ation lies between those of a maser and a laser, we hereby 
declare the device to be a Josephson STAR-emitter, for 
stimulated terahertz amplified radiation emitter. 
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